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Abstract 

Density  functional  theory  has  been  applied  to  a  study  of  the  electronic  structure  of  the  ideally-terminated,  relaxed 
and  H-saturated  (0001)  surfaces  of  P-Si3N4  and  to  that  of  the  bulk  material.  For  the  bulk,  the  lattice  constants  and 
atom  positions  and  the  valence  band  density  of  states  are  all  in  good  agreement  with  experimental  results.  A  band 
gap  of  6.7  eV  is  found  which  is  in  fair  accord  with  the  experimental  value  of  5. 1-5.3  eV  for  H-free  Si3N4.  Using  a 
two-dimensionally-periodic  slab  model,  a  n-bonding  interaction  is  found  between  threefold-coordinated  Si  and  two- 
fold-coordinated  N  atoms  in  the  surface  plane  leading  to  n  and  n*  surface-state  bands  in  the  gap.  A  surface-state  band 
derived  from  s-orbitals  is  also  found  in  the  gap  between  the  upper  and  lower  parts  of  the  valence  band.  Relaxation 
results  in  displacements  of  surface  and  first-underlayer  atoms  and  to  a  stronger  7i-bonding  interaction  which  increases 
the  7t-7t*  gap.  The  relaxed  surface  shows  no  occupied  surface  states  above  the  valence  band  maximum,  in  agreement 
with  recent  photoemission  data  for  a  thin  Si3N4  film.  The  n*  band,  however,  remains  well  below  the  conduction  band 
minimum  (but  well  above  the  Fermi  level).  Adsorbing  FI  at  all  dangling-bond  sites  on  the  ideally-terminated  surface  and 
then  relaxing  the  surface  and  first  underlayer  leads  to  smaller,  but  still  finite,  displacements  in  comparison  to  the  clean 
relaxed  surface.  This  surface  is  more  stable,  by  about  3.67  eV  per  H,  than  the  clean  relaxed  surface. 
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1.  Introduction 

Silicon  nitride  fSi3N4)  is  an  important  electronic 
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coating  for  magnetic  storage  devices.  In  contrast 
to  the  situation  for  other  such  materials,  the  sur¬ 
face  electronic  structure  of  Si3N4  is  not  well  under¬ 
stood  and  appears  to  depend  on  the  details  of  how 
the  nitride  is  formed.  For  example,  ultraviolet  pho¬ 
toemission  spectroscopy  (UPS)  data  for  a  Si(l  11)- 
(7  x  7)  surface  nitrided  by  reaction  with  NFI3  at 
high  temperature  [1]  show  a  surface  state  at 
1.1  eV  above  the  valence  band  maximum  (VBM). 
On  the  other  hand,  a  disordered  (polycrystalline) 
Si3N4  film,  grown  on  Si(  100)  by  low-pressure 
chemical  vapor  deposition  and  chemically  cleaned 
in  situ  in  ultrahigh  vacuum  (UHV),  shows  no  evi¬ 
dence  of  surface  states  near  the  VBM  [2].  As  noted 
elsewhere  [2]  the  vibrational  spectrum  of  a  Si  sur¬ 
face  nitrided  by  reaction  with  NH3  is  also  mark¬ 
edly  different  from  that  of  bulk  Si3N4.  Charge 
trapping  at  Si3N4  surfaces,  which  may  be  defect-re¬ 
lated,  is  a  further  important  but  not  well-under¬ 
stood  aspect  of  this  material1  [Refs.  [3-6]  and 
works  cited].  One  reason  for  the  lack  of  informa¬ 
tion  about  atomically-clean  and  stoichiometric 
surfaces  of  bulk  Si3N4  is  the  fact  that  a  procedure 
for  preparing  such  surfaces  in  UHV  has  only  re¬ 
cently  been  reported  [2]. 

The  bulk  electronic  structure  of  Si3N4  has  been 
studied  extensively  [7-17]  using  theoretical  meth¬ 
ods,  and  many  groups  have  performed  computa¬ 
tional  studies  of  either  the  chemisorption  of 
atomic  N  [18-20]  or  the  initial  decomposition  of 
NH3  (e.g.,  Refs.  [21-27])  on  Si  surfaces.  However, 
less  such  work  has  been  reported  for  the  surfaces 
of  bulk  Si3N4.  Bagatur’yants  et  al.  [28]  performed 
an  ab  initio  (RHF/MP2)  cluster  study  of  chemical 
reactions  that  occur  on  the  p-Si3N4(000 1)  surface 
during  growth  by  chemical  vapor  deposition.  In 
particular,  relaxation  of  the  clean  surface  was 
found  to  produce  >Si=N-  structures  as  a  means 
of  eliminating  dangling  orbitals  on  the  ideally-ter¬ 
minated  surface  (see  below).  Zhao  and  Bachlech- 
ner  [11]  considered  a-  and  |3-Si3N4  and  found 
that  Si-N  bonds  are  more  covalent  on  the  unre¬ 
laxed  surfaces  than  in  the  bulk. 


1  The  recurring  reference  to  “SfiNLt”,  instead  of  in 

Ref.  [6]  is  a  misprint. 


In  this  work  we  will  examine  the  physical  and 
electronic  structure  of  the  ideally-terminated,  the 
relaxed  and  the  H-saturated  P-Si3N4(000 1) 
surface  using  a  two-dimensionally-periodic  slab 
model  in  an  ab  initio  density  functional  theory 
(DFT)  calculation.  This  will  allow  the  surface 
band  structure  to  be  obtained,  which  is  not  pro¬ 
vided  in  a  treatment  of  an  isolated  cluster.  The  re¬ 
sults  will  also  provide  an  understanding  of  the 
nature  of  the  P-Si3N4(000 1)  surface  states  which 
have  not  been  addressed  in  detail  previously. 

2.  Computational  details 

All  calculations,  except  where  noted,  were  per¬ 
formed  using  the  CRYSTAL2003  suite  of  pro¬ 
grams  [29,30]  which  uses  Bloch  functions 
constructed  from  localized  Gaussian  basis  sets. 
The  Si  basis  set  was  the  65-11 1G*  set  developed 
by  D’Arco  et  al.  [31]  for  partially-ionic  Si  com¬ 
pounds,  with  the  exponents  of  the  outer  sp-  and 
d-shells  re-optimized  for  Si3N4.  The  results  were 
1.249,  0.5477  and  0.2834  for  the  sp-shells  and 
0.5437  for  the  d-shell  (all  in  units  of  bohr-2). 
The  N  basis  set  was  the  same  7-31 1G*  set  used 
for  GaN  [32],  with  exponents  of  the  outer  sp- 
and  d-shells  re-optimized  for  Si3N4.  The  results 
were  0.4302  and  0.1357  for  the  sp-shells  and 
0.5976  for  the  d-shell. 

Optimization  of  the  basis  sets  and  the  unit-cell 
lattice  constants  was  done  using  the  LoptCG 
script2  which  computes  gradients  numerically 
using  CRYSTAL98,  followed  by  the  application 
of  a  conjugate-gradient  algorithm  to  locate  the  to¬ 
tal-energy  minimum.  Optimization  of  the  positions 
of  individual  atoms  in  the  unit  cell  or  on  the  sur¬ 
face  were  done  using  CRYSTAL2003  which  ob¬ 
tains  energy  gradients  analytically.  The  DFT 
calculations  used  the  hybrid  B3LYP  functional 
which  has  been  shown  [33,34]  to  provide  more 


2  The  LoptCG  script  was  obtained  at  http://www.crystal. 
unito.it/.  At  present,  the  publicly-available  version  of  the  script 
works  only  with  CRYSTAL98.  The  latest  version,  CRYS- 
TAL2003,  can  generate  gradients  analytically  but  only  for  the 
positions  of  individual  atoms,  not  for  basis-set  parameters  or 
lattice  constants. 
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accurate  band  gaps  than  do  purely  ab  initio  func¬ 
tionals.  The  exchange-correlation  potential  was 
determined  by  numerical  integration  [29]  rather 
than  by  fitting  with  Gaussian  functions  [35]  as 
was  done  in  CRYSTAL98.  Due  to  the  presence 
of  singly-occupied  dangling  orbitals  on  the  ide¬ 
ally-terminated  (0001)  surface  (see  below),  spin- 
unrestricted  DFT  was  used  in  all  calculations 
involving  clean  surfaces.  Si3p  (N2p)  electrons  on 
surface  sites  with  dangling  orbitals  were  arbitrarily 
assigned  initial  spins  of  +1/2  (—1/2)  to  start  the 
first  self-consistent  field  (SCF)  cycle.  A  shrinking 
factor  of  six  was  used  for  all  reciprocal-lattice  vec¬ 
tors  in  the  k-point  interpolations. 


3.  Results  and  discussion 

3.1.  Bulk  properties 

The  physical  and  electronic  structure  of  the 
bulk  were  analyzed  as  a  necessary  preliminary  to 
addressing  surface-related  issues.  For  bulk  (3- 
Si3N4  the  P63/m  structure  (space  group  no.  176) 
was  used  [8,9],  as  shown  in  Fig.  1.  The  lattice  con¬ 
stants  and  atomic  positions  were  optimized  as  de¬ 
scribed  above,  and  the  results  (Table  1)  are  found 
to  be  in  good  accord  with  experimental  values  [36] 
and  with  previous  theoretical  values. 

Fig.  2  shows  the  computed  VB  density  of  states 
(DOS)  together  with  results  from  UPS  [2]  and 
from  X-ray  photoemission  spectroscopy  (XPS) 
[37,38]  for  H-free  Si3N4.  The  orbital  composition 


Fig.  1.  (a)  Unit  cell  for  the  bulk  P63/m  lattice  structure  of  f- 
Si3N4,  with  atoms  labelled  as  in  Ref.  [28],  (b)  The  ideal 
(unrelaxed)  surface  viewed  along  the  surface  normal.  Solid 
(open)  symbols  refer  to  surface  (first-underlayer)  atoms,  and  the 
various  inequivalent  atom  types  are  labelled.  Bonds  are  shown 
between  surface  atoms,  with  double  lines  indicating  the  ir-bonds 
that  form  when  relaxation  occurs.  For  clarity,  bonds  involving 
underlayer  atoms  are  not  shown.  The  unit-cell  edges  are  shown 
in  both  panels. 

of  the  various  parts  of  the  VB,  which  will  be 
important  in  the  following  discussion  of  surface 
states,  has  been  analyzed  previously  [14-16],  The 
lower  part  of  the  VB  (from  about  —19  to  —22  eV 
in  Fig.  2)  consists  mainly  of  N  s-orbitals  with  a 
small  contribution  from  Si  s-orbitals.  In  the  upper 
part  of  the  VB,  states  near  the  VBM  arise  mainly 
from  non-bonding  N  p-orbitals  oriented  perpen¬ 
dicular  to  the  plane  of  the  N-Si-N  bond,  while 


Table  1 

Observed  and  calculated  lattice  parameters  for  p-Si3N4  (P63/m  structure) 


Reference 

Lattice  constant  (A) 

Si(6h)a 

N2(6h)a 

a 

c 

X 

y 

X 

y 

This  workb 

7.655 

2.918 

0.1738 

0.7679 

0.3303 

0.0299 

Ref.  [7]b 

7.6473 

2.9171 

0.1751 

0.7689 

0.3302 

0.0307 

Ref.  [9]b 

7.562 

2.893 

0.1741 

0.7677 

0.3302 

0.0298 

Ref.  [10]b 

7.622 

2.910 

Ref.  [15]b 

7.61 

2.91 

Ref.  [36]c 

7.608 

2.911 

0.1733 

0.7694 

0.3323 

0.0314 

“  The  Nl(2c)  coordinates  were  fixed  at  x  =  1/3,  y  =  2/3.  The  --coordinate  for  all  atoms  was  fixed  at  z  =  1/4.  Flere  “x”  and  “y”  are 
fractions  of  the  lattice  constant  “a”,  and  “z”  is  a  fraction  of  “c”. 
b  Theoretical  results.  See  references  for  details  of  the  calculations. 
c  Experimental  results. 
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Binding  Energy  (eV) 

Fig.  2.  Comparison  of  photoemission  data  with  the  calculated 
VB  DOS  for  bulk  p-Si3N4.  (a)  Hell  UPS  (hv  =  40.8  eV)  from 
Ref.  [2],  (b)  MgKoc  XPS  (hv  =  1253.6  eV)  from  Ref.  [38],  (c) 
Calculated  VB  DOS  (this  work).  For  all  plots,  binding  energy  is 
referenced  to  the  calculated  Fermi-level  position,  which  falls 
near  mid-gap.  The  UPS  and  XPS  data  were  shifted  in  energy  to 
align  the  peak  near  —  1 3  eV  with  the  corresponding  feature  in 
the  calculated  DOS. 


deeper-lying  states  are  due  to  Si-N  a-bonding 
orbitals.  The  bottom  of  the  upper  part  of  the  VB 
(near  — 13  eV  in  Fig.  2)  is  derived  largely  from  Si 
s-orbitals. 

The  width  and  overall  shape  of  the  DOS  agree 
well  with  experiment,  as  do  the  results  of  previous 
calculations.  The  different  appearance  of  the  upper 
part  of  the  VB  in  UPS  and  in  XPS  results  from  the 
increase  in  the  relative  photoionization  cross-sec¬ 
tion  for  s-  vs.  p-orbitals  at  higher  incident  photon 
energies  [39].  One  notes  in  particular  the  close  sim¬ 
ilarity  between  the  UPS  results,  which  were  ob¬ 
tained  for  an  atomically-clean  polycrystalline 
surface  under  conditions  of  high  surface  sensitivity 
[2],  and  the  corresponding  bulk  DOS  of  crystalline 


(3-Si3N4.  This  supports  the  previous  conclusion 
regarding  the  absence  of  surface  states  above  the 
VBM  for  the  Si3N4  samples  discussed  in  Ref.  [2]. 
Finally,  a  band  gap  of  6.7  eV  was  found  in  the 
present  calculation  vs.  the  experimental  value  of 
5.1-5.3eV  reported  [40  42]  for  FI-free  Si3N4.  As 
noted  elsewhere  [33],  use  of  the  B3LYP  functional 
typically  leads  to  a  band  gap  that  is  somewhat  lar¬ 
ger  than  the  experimental  result. 

3.2.  The  clean  (0001)  surface 

3.2.1.  The  ideally-terminated  surface 

The  surface  was  modeled  using  a  two-dimen- 
sionally  periodic  slab  which  was  six  layers  thick. 
Six  layers  are  equivalent  to  three  unit  cells  (cf. 
Fig.  1)  stacked  along  the  c-axis.  This  is  sufficiently 
thick  that  no  surface-related  features  appeared  in 
the  VB  partial  density  of  states  (PDOS)  for  the 
middle  of  the  slab  and  that  the  Mulliken  charges 
on  atoms  in  the  middle  of  the  slab  were  virtually 
the  same  as  for  the  bulk  material.  Thus,  six  layers 
was  considered  to  be  sufficient  to  separate  the  two 
slab  surfaces  from  each  other.  Since  the  present 
work  uses  localized,  rather  than  plane-wave,  basis 
sets  it  is  not  necessary  to  have  periodically- 
repeated  slabs  along  the  surface  normal. 

Referring  to  Fig.  1,  three  Si  and  four  N  atoms 
lie  in  the  ideally-terminated  surface  plane.  Three 
of  the  N’s  (“terminal”,  labelled  “Nt”  as  in  Ref. 
[28])  are  bonded  to  only  two  Si’s,  one  in  the  surface 
plane  and  one  in  the  first  underlayer.  Therefore, 
each  Nt  has  a  dangling  orbital  (in  addition  to  the 
doubly-occupied  non-bonding  lone-pair  orbital). 
The  fourth  N  (“central”,  labelled  “Nc”)  is  bonded 
to  three  Si’s  in  the  surface  plane  and  thus  has  no 
singly-occupied  dangling  orbital.  All  three  surface 
Si’s  are  bonded  to  only  three  N’s,  and  each  has  one 
dangling  orbital.  In  the  present  calculations,  the 
dangling  bonds  on  the  identical  “back”  surface 
of  the  slab  were  all  terminated  with  H  atoms,  for 
which  the  standard  Pople  6-3 1G**  basis  set  was 
used. 

In  qualitative  agreement  with  previous  results 
[11],  the  Si-N,  bond  on  the  ideally-terminated  sur¬ 
face  is  found  to  be  somewhat  less  ionic  than  in  the 
bulk.  The  present  results  give  +1.99  (—1.47)  for  the 
partial  charge  on  a  surface  Si(Nt)  atom  vs.  +2.15 
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(—1.63)  for  the  corresponding  bulk  values.  On  the 
other  hand,  and  also  in  accord  with  Ref.  [11],  the 
partial  charge  on  a  surface  Nc  atom  (—1.59)  is  al¬ 
most  identical  to  that  of  a  bulk  Nc  (—1.56).  The 
Mulliken  population  analysis  also  indicates  a  sig¬ 
nificant  n-bonding  interaction  between  surface  Si 
and  Nt  sites,  even  in  the  absence  of  relaxation, 
which  affects  the  structure  of  surface  states  on 
the  ideally-terminated  surface  (see  below).  This 
interaction  leads  to  a  Si-Nt  overlap  population 
of  0.392  at  the  surface  vs.  a  value  of  0.271  in  the 
bulk. 

Figs.  3a  and  4a  show  band-structure  results 
which  indicate  the  presence  of  surface  states  in 
the  lower  part  of  the  VB  as  well  as  near  the  bulk 
VBM.  The  bulk  VBM  was  determined  by  identify¬ 
ing  the  highest-energy  band  that  is  composed  lar¬ 
gely  of  orbitals  centered  on  atoms  not  in  the 
surface  plane.  For  clarity,  only  states  within  a 
few  eV  of  the  band  edges  are  included  in  Fig.  3. 
Examination  of  the  DOS  (not  shown)  indicates 
that  the  occupied  surface  states  are  derived  mainly 
from  the  surface  Nt,  with  a  smaller  contribution 


from  surface  Si,  and  that  the  occupancies  are 
somewhat  different  for  “up”  and  “down”  spins 
as  expected  for  states  with  singly-occupied,  dan¬ 
gling-bond  character.  The  occupied  states  near 
the  VBM  have  a  7r-bonding  character.  The  corre¬ 
sponding  7t*-antibonding  states  are  seen  at 
approximately  mid-gap  in  Fig.  3a. 

3.2.2.  The  relaxed  (0001)  surface 

The  cluster  results  of  Bagaturyants  et  al.  [28] 
indicate  that  surface  relaxation  takes  the  form  of 
7t-bond  formation  between  Si-Nt  nearest-neigh- 
bors  which  eliminates  singly-occupied  dangling 
orbitals.  In  the  present  work,  relaxation  was  trea¬ 
ted  by  starting  with  the  ideally-terminated  surface 
and  allowing  atoms  in  the  surface  layer  and  the 
first  underlayer  to  move  freely  during  optimiza¬ 
tion.  Relaxation  was  found  to  lower  the  total  en¬ 
ergy  by  3.12  eV  per  (lxl)  surface  unit  cell. 

Table  2  summarizes  the  structural  results  which 
are  in  qualitative  agreement  with  those  of  previous 
work  [28].  Relaxation  results  in  an  inward  (out¬ 
ward)  displacement  of  surface  Si(Nc)  atoms  and 


(a)  (b) 


K  r  MKK  r  MK 


Fig.  3.  Partial  band  structure  diagram  for  (a)  the  ideally-terminated  and  (b)  the  relaxed  p-Si3N4(000 1)  surface  for  k-vectors  in  the 
plane  of  the  two-dimensionally-periodic  six-layer  slab.  The  zero  of  energy  is  the  Fermi  level.  States  in  the  VB  below  about  —5.5  eV  and 
in  the  CB  above  about  6  eV  are  omitted  (cf.  Fig.  2).  The  arrows  mark  the  approximate  position  of  the  "bulk”  VBM  (see  text).  The  n 
and  n*  surface-state  bands  involving  surface  Si  and  N,  sites  are  labeled  in  panel  (a).  For  the  it*  bands  of  panel  (a),  solid  and  dashed 
lines  (barely  distinguishable  on  the  energy  scale  of  the  figure)  show  the  “spin-up”  and  “spin-down”  bands. 
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(a) 


(b) 


Fig.  4.  Similar  to  Fig.  3  but  showing  the  gap  between  the  lower  and  upper  parts  of  the  VB  (cf.  Fig.  2).  Valence  states  lying  above 
— 12  eV  have  been  omitted  for  clarity. 


Table  2 


Structural  parameters  computed  for  the  relaxed  surface3 


Si-Nc 

Si=Nt 

Si-Nt 

Nt-Si 

Sz(Si) 

Sz(Nc) 

5r(Nt) 

Surface 

1.758 

1.601 

1.687b 

1.71 1° 

—0.16 

+0.27 

+0.02 

First  underlayer*1 

1.738 

1.774 

+0.08 

+0.04 

+0.02 

Bulk 

1.743 

1.748 

Cluster  [28f 

1.75 

1.61 

-0.24 

+0.07 

+0.10 

a  All  distances  are  in  Angstroms.  Si-Nc,  etc.  are  bond  lengths.  The  8r  quantities  are  displacements  along  the  surface  normal  relative 
to  the  atomic  positions  on  the  ideally-terminated  surface.  Omitted  values  are  either  not  applicable  or  are  zero  by  definition. 
b  This  is  the  back-bond  distance  between  a  surface  Si  and  a  first- underlayer  N,  (cf.  Fig.  1). 
c  This  is  the  back-bond  distance  between  a  surface  N,  and  a  first-underlayer  Si  (cf.  Fig.  1). 
d  These  distances  refer  to  bonds  within  the  first  underlayer. 
e  Results  for  the  3-DAS  cluster  model  discussed  in  Ref.  [28], 

a  slight  outward  displacement  of  the  Nt,  all  of 
which  are  confined  mainly  to  the  surface  layer. 

This  results  in  the  Nc  atom  lying  about  0.43  A 
above  the  plane  of  the  three  surface  Si  atoms  to 
which  it  is  bonded.  In-plane  (lateral)  displacements 
of  surface  atoms  occur  mainly  for  the  Si’s,  with  a 
magnitude  of  about  0.22  A.  Those  for  the  N  atoms 
are  small  (^0.04  A  for  Nt  and  for  Nc).  In  con¬ 
trast  to  the  atomic  positions,  the  atomic  charges 
on  the  relaxed  surface  (+2.00,  —1.41  and  —1.49 
for  Si,  Nt  and  Nc,  respectively)  are  nearly  the  same 
as  for  the  ideally-terminated  surface  (given  above). 


Formation  of  the  Si=Nt  Tt-bond  is  recognizable  in 
the  reduced  interatomic  distance  (1.601  A)  and 
increased  Mulliken  overlap  population  (0.475)  vs. 
the  bulk  Si-Nt  o-bond  (1.748  A  and  0.271,  respec¬ 
tively).  The  back-bonds  from  surface  to  underlayer 
atoms  are  also  somewhat  reduced  in  length  and 
increased  in  overlap  population  in  comparison  to 
bulk  bonds.  On  the  other  hand,  the  Si-Nt  cr-bond 
within  the  first  underlayer  is  slightly  weakened 
(i.e.,  lengthened). 

The  surface  core-level  shift  (SCLS)  in  binding 
energy  (BE)  for  the  N  Is  and  Si2p  levels  at  differ- 
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Table  3 

Relative  binding  energy  shifts  calculated  for  Ihe  N  Is  and  Si2p  core  levels  for  different  forms  of  the  SuN^OOO  1)  surface3 


Relaxed 

Ideally-term. 

H-saturated 

Nls 

Si2p 

Nls 

Si2p 

Nls 

Si2p 

Nt  bulk 

0 

0 

0 

Nc  bulkb 

-0.47 

-0.44 

-0.34 

N,  surface 

-0.87 

-0.73 

+0.26 

Nc  surface 

Si  bulk 

+  1.68 

0 

+  1.31 

0 

-0.14 

0 

Si  surface3 

+0.68 

+0.52 

-0.27 

a  All  quantities  are  in  eV,  relative  to  either  N,  bulk  or  Si2p  bulk.  Negative  (positive)  values  indicate  a  decrease  (increase)  in  BE. 
b  The  slightly  different  Nc  bulk  values  reflect  the  level  of  precision  in  determining  the  BE  shifts. 
c  A  ~0.2  eV  crystal  field  splitting  of  the  clean-surface  Si2p  is  neglected. 


ent  sites  was  estimated  for  comparison  with  X-ray 
photoemission  spectroscopy  (XPS)  results  [2].  This 
was  done  simply  by  using  the  relative  orbital  ener¬ 
gies  for  different  inequivalent  atomic  species.  Such 
an  “initial-state”  (Koopman’s  theorem)  approach 
neglects  differences  in  core-hole  screening  [43]  but 
should  yield  reasonable  results  for  SCLSs  (as  op¬ 
posed  to  absolute  BEs).  Table  3  summarizes  the  re¬ 
sults.  The  Si  and  N,  SCLS  values  exhibit  the  trend 
expected  [43]  for  partially-ionic  materials  in  which 
cation  (anion)  core  levels  shift  to  higher  (lower)  BE 
due  to  the  decrease  in  the  number  of  oppositely- 
charged  nearest-neighbor  sites  at  the  surface  vs. 
in  the  bulk.  The  Nt  SCLS  is  too  small  to  be  de¬ 
tected  in  the  data  presented  in  Ref.  [2]  given  the 
limited  instrumental  resolution  and  surface  sensi¬ 
tivity  and  the  N  Is  line  width.  It  is  noted  that 
core-level  line  widths  in  partially-ionic  materials 
such  as  Si3N4  receive  a  significant  contribution 
from  zero-point  lattice  vibration  [44]  in  addition 
to  the  effects  of  any  sample  imperfections.  Like¬ 
wise,  the  estimated  +0.68  eV  shift  for  the  surface 
Si2p  was  not  resolved  [2]  in  the  Zr  Mc-cxcited 
(hv  =  151.6  eV)  Si2p  XPS  data.  In  this  case,  the 
instrumental  resolution  and  surface  sensitivity 
were  much  better  than  in  the  N  Is  data,  and  a  shift 
of  this  magnitude  should  have  been  observable. 
The  results  suggest  that  the  magnitude  of  the 
SCLS  is  overestimated  in  Table  3,  at  least  for  the 
Si2p. 

A  substantial  shift  to  higher  BE  is  found  for  the 
Ncls  level  at  the  relaxed  surface.  Relative  to  the 
bulk  Nt,  which  contributes  most  of  the  N 1  s  inten¬ 
sity  for  MgKa  excitation  (hv  -  1253.6  eV),  the  sur¬ 


face  Ncls  is  shifted  1.68  eV  to  higher  BE.  This 
agrees  well  (perhaps  fortuitously  so)  with  a  sur¬ 
face-related  satellite  observed  [2]  at  1.7  eV  below 
the  main  N  Is  XPS  peak  for  clean  Si3N4  surfaces.1 * 3 
Each  Nc  on  the  relaxed  surface  is  in  a  configura¬ 
tion  that  can  be  depicted  (cf.  Fig.  lb)  as  Nc(- 
Si=Nt-)3,  and  the  large  positive  SCLS  appears  to 
be  correlated  with  the  strength  of  the  Si=Nt  n- 
bond.  The  positive  SCLS  is  greatest  for  the  relaxed 
surface,  smaller  for  the  ideally-terminated  surface 
and  absent  (i.e.,  a  small  negative  SCLS)  on  the 
H-saturated  surface  (see  below). 

Figs.  3b  and  4b  show  the  band  structure  dia¬ 
gram  for  the  relaxed  surface.  The  surface  states 
evident  in  Figs.  3a  and  4a  for  ideal  termination 
are  modified  by  the  relaxation,  which  increases 
the  Si-Nt  interaction  and  shifts  the  ir-bonding 
states  down  in  energy  and  the  7i*-antibonding 
states  farther  up  into  the  gap.  The  empty  bands 
in  Fig.  3  represent  virtual  states  and  cannot  be  re¬ 
lated  quantitatively  to  observed  inter-band  transi¬ 
tions.  However,  the  results  suggest  that  transitions 
at  energies  less  than  the  bulk  band  gap  should  be 
observable  in  surface-sensitive  electron  energy  loss 
spectroscopy.  Such  an  excitation  has  been  re¬ 
ported  [2]  at  ~4.6  eV  which  might  be  related  to 


1  Recent  work  [45]  has  found  an  Nls  SCLS  of  0.7  eV  to 

higher  BE  for  thin  films  prepared  by  plasma  nitridation  of  a 

Si(  111)  surface.  Since  different  methods  were  used  to  prepare 
the  surfaces  studied  in  Refs.  [2,45],  it  is  unknown  whether  the 
physical  origin  of  the  satellite  is  the  same  in  both  cases. 
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Table  4 

Structural  parameters  computed  for  the  H-saturated  surface3 


HSi-N0 

HSi-N,H 

HSi-N, 

HN,-Si 

8z(SiH) 

6z(Nc) 

8z(NtH) 

Surface 

1.744 

1.726 

1.749b 

1.730° 

+0.02 

+0.17 

-0.04 

First  underlayer*1 

1.751 

1.745 

+0.02 

+0.12 

0.00 

Bulk 

1.743 

1.748 

Cluster  [28]e 

1.76 

1.74 

-0.01 

+0.07 

+0.04 

a  See  Table  2,  footnote  (a).  In  the  present  table,  the  optimized  Si-H  and  Nt-H  distances  are  1.469  and  1.013  A,  respectively. 
b  This  is  the  back-bond  distance  between  a  surface  SiH  and  a  first-underlayer  Nt. 

0  This  is  the  back-bond  distance  between  a  surface  NtH  and  a  first-underlayer  Si. 
d  These  distances  refer  to  bonds  within  the  first  underlayer.  No  H  is  present  in  this  layer. 
e  Results  for  the  (3-DAS)H6  model  discussed  in  Ref.  [28], 


the  7i  — >  7i*  transition  of  the  Si=N,  bond.  In  the  re¬ 
gion  between  the  upper  and  lower  parts  of  the  VB 
(Fig.  4),  relaxation  shifts  the  surface  states  to  a  po¬ 
sition  a  few  tenths  of  an  eV  above  the  bulk  maxi¬ 
mum  of  the  lower  part  of  the  VB. 

3.3.  The  H-saturated  (0001)  surface 

Next  we  consider  the  surface  that  results  when 
the  singly-occupied  dangling  bonds  (three  on  Si 
atoms  and  three  on  Nt  sites)  on  the  ideally-termi¬ 
nated  surface  are  “tied  up”  with  H.  In  the  geo¬ 
metry  optimization,  the  positions  of  the  six  H 
adatoms  and  of  the  surface  and  first-underlayer 
atoms  were  allowed  to  vary  freely.  As  in  the  case 
of  the  other  calculations  described  above,  the 
“back”  surface  of  the  slab  was  terminated  by  H 
atoms  in  fixed  positions.  The  structural  results 
given  in  Table  4  show  that  the  interatomic  dis¬ 
tances  for  the  H-saturated  surface  are  much  closer 
to  the  corresponding  bulk  values,  and  the  8z  values 
generally  smaller,  than  was  the  case  for  the  clean 
relaxed  surface.  As  expected,  the  SCLS  values  for 
the  H-saturated  surface  (Table  3)  are  much  less 
than  for  either  clean  surface. 

The  energy  difference,  AEH,  between  the  (clean) 
relaxed  and  H-saturated  surfaces  is  given  by 

£(relax.)  +  6£(H)  =  £(H-sat.)  +  3£Zp(Si-H) 

+  3iizp(N-H)  +  A £fi 

where  iffrelax.)  and  A(H-sat.)  are  the  geometry- 
optimized  total  energies  of  the  respective  surfaces 
and  E( H)  is  the  energy  of  the  free  H  atom.  In  this 
definition,  a  positive  A£H  indicates  an  exothermic 


process.  The  calculation  uses  H  atoms,  rather  than 
molecular  H2,  since  it  is  standard  practice  for  H 
chemisorption  experiments  on  non-metals  to  be 
performed  using  a  hot  (~2000  K)  tungsten  fila¬ 
ment  to  dissociate  the  H2  molecules.  The  E/y> 
terms  are  the  zero-point  vibrational  energies  of 
the  Si-H  and  N-H  bonds.  Any  changes  in  the  nor¬ 
mal  modes  of  the  slab  itself  that  may  result  from 
adsorption  are  ignored  in  applying  the  zero- 
point  corrections.  For  the  basis  set  used  here, 
F(H)  = —13.558  eV  (close  to  the  experimental 
value  of  -13.606  eV),  and  FzpISi-H)  =  0.130  eV 
and  £zp(N-H)  =  0.21 1  eV  based  on  spectroscopic 
data  [2]  for  the  bond-stretching  normal  modes. 
The  final  result  is  A/:H  =  +22.021  eV  per  surface 
unit  cell,  or  an  average  of  3.67  eV  per  adsorbed 
H,  which  indicates  a  substantial  stabilization  vs. 
the  clean  relaxed  surface.  This  can  be  compared 
with  typical  molecular  Si-H  and  N-H  bond  ener¬ 
gies  of  about  3.0  and  4.0  eV,  respectively  [46], 
The  stabilization  probably  derives  in  part  from 
the  reduction  in  strain  that  accompanies  elimina¬ 
tion  of  the  Si=Nt  7t-bond,  in  which  the  Si  is  forced 
into  an  approximately  planar,  threefold-coordi¬ 
nated  configuration. 


4.  Summary  and  conclusions 

Density  functional  theory,  with  the  B3LYP 
functional  and  all-electron  Gaussian  basis  sets, 
has  been  applied  to  a  determination  of  the  elec¬ 
tronic  structure  of  the  ideally-terminated,  relaxed 
and  H-saturated  (0001)  surfaces  of  p-Si3N4  as 
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well  as  to  that  of  the  bulk  material.  A  two-dimen- 
sionally-periodic  slab  model  was  used  for  the  sur¬ 
face  calculations.  The  results  are  as  follows. 

(1)  For  the  bulk,  the  lattice  constants  and  atom 
positions  and  the  VB  DOS  are  all  in  good 
agreement  with  experimental  results.  A  band 
gap  of  6.7  eV  is  found  which  is  reasonably 
close  to  the  experimental  value  of  5.1- 
5.3  eV  for  H-free  Si3N4. 

(2)  For  the  ideally-terminated  surface  a  Ji-bond- 
ing  interaction  is  found  between  threefold- 
coordinated  Si  and  twofold-coordinated  N 
atoms  in  the  surface  plane  leading  to  ji  and 
n*  surface-state  bands  in  the  gap.  A  sur¬ 
face-state  band  derived  from  s-orbitals  is  also 
found  in  the  gap  between  the  upper  and 
lower  parts  of  the  VB. 

(3)  Relaxation  leads  to  displacements  of  surface 
and  first-underlayer  atoms  and  to  a  stronger 
7t-bonding  interaction  which  increases  the  7t- 
7i*  gap.  The  relaxed  surface  shows  no  occu¬ 
pied  surface  states  above  the  VBM  (in  agree¬ 
ment  with  recent  UPS  data  [2]  for  a  thin 
Si3N4  film).  The  tt*  band,  however,  remains 
well  below  the  CBM  (but  well  above  the 
Fermi  level)  and  might  be  the  final  state  in 
a  sub-band-gap  excitation  observed  [2]  in 
ELS. 

(4)  A  high-BE  N  Is  satellite  seen  in  XPS  is  iden¬ 
tified  with  the  threefold-coordinated  (Nc)  site 
on  the  relaxed  surface  on  the  basis  of  a 
Koopman’s- theorem  estimate  of  the  SCLS. 
This  corresponds  to  an  experimentally- 
observed  Nls  satellite  [2],  Smaller  SCLSs 
are  predicted  for  Si2p  and  Ntls  and  have 
not  yet  been  detected  experimentally. 

(5)  Adsorbing  H  at  all  dangling-bond  sites  on 
the  ideally-terminated  surface  and  then 
relaxing  the  surface  and  first  underlayer  leads 
to  smaller,  but  still  finite,  displacements  in 
comparison  to  the  bare  surface.  The  hydro¬ 
genated  surface  is  more  stable,  by  about 
3.67  eV  per  H,  than  the  clean  relaxed  surface. 

These  results  provide  a  structural  model  for  the 
Si3N4  surface  which  will  be  useful  as  a  starting 
point  in  interpreting  data  from  chemisorption 


experiments  and  in  understanding  the  properties 
of  defects  on  this  surface. 
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